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Abstract 

We make a critical reappraisal of 'unbounded-from-below' (UFB) constraints in the 
MSSM and i?-parity violating models. We explain why the 'traditional' UFB bounds 
are neither necessary nor sufficient and propose, instead, a sufficient condition which 
ensures that there are no local minima along the flat directions. This conservative 
(but meaningful) condition divides the parameter space into regions which are al- 
lowed, regardless of cosmology, and regions in which cosmology is expected to play 
a major role. We study both conditions at low tan/3 and obtain analytic approxi- 
mations to the UFB bounds for all low tan (3 (< 15). Finally we show that i?-parity 
violation just below current experimental limits avoids UFB problems by lifting the 
dangerous flat directions. 
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1 Introduction 



Metastability is an important issue for supersymmetric models because of the existence 
of directions which are flat in the limit of unbroken supersymmetry. Supersymmetry 
breaking can result in the physical vacuum being metastable; an uncomfortable situation 
which some authors believe is grounds for a model to be rejected. Others argue that 
a model which is metastable is fine, provided that the vacuum could not have decayed 
within the lifetime of the universe (see Refs. |IJ]-[||] for examples of both). This particular 
question is unlikely to be resolved until more is known about early cosmology. In the 
meantime, it is important to ask whether metastability is actually the best criterion for 
deciding if a flat direction is going to be troublesome for cosmology. In fact in this paper 
we will argue that it is not. We shall also be considering a non-cosmo logical solution to 
the problems posed by flat directions - i?-parity violation. 

To address these issues we shall re-examine F and D-flat directions in the Minimal 
Supersymmetric Standard Model (MSSM) and in .R-parity violating models. D-flat di- 
rections can be driven negative by the soft supersymmetry breaking terms leading to 
undesirable charge and colour breaking minima |], [| 0, || . When the direction is also F- 
flat, particularly strong constraints on parameter space arise which go by the misnomer 
of 'Unbounded From Below' (UFB) constraints J| [3], Provided that all the super- 
symmetry breaking mass-squared terms are positive at the high (we shall assume GUT) 
scale, what actually happens is that a minimum forms radiatively at a scale of order 
few x mw/hb- Imposing the condition that this minimum is not the global one results in 
what we shall refer to as the 'traditional' UFB bound. 

Following the original work on charge and colour breaking vacua [jl], |2|], such bounds 
have been analysed directly using numerical runnings of the MSSM [|3], [|]. Although it 
gives good accuracy, this type of analysis is less than transparent. Instead, we will work 
with analytic solutions to the one-loop RGEs in the approximation that the bottom and 
tau contributions are negligible (valid for all low tan/3). These solutions will enable us 
to derive traditional UFB limits (on m ) for all low tan/3 (as a function of M1/2, tan/3 
and and A ). Thus by sacrificing some accuracy we gain a lot of flexibility and also 
some insight into how the potentials along the flat directions behave as we change the 
parameters. (We should stress that we are not aiming to improve on the accuracy of 
previous calculations, indeed our estimates will be considerably worse, being accurate to 
maybe 15% at best.) Analytic solutions also allow us to address, in more detail than is 
usually possible, two questions: 

• Are the traditional UFB bounds either necessary or sufficient? 

• Can R-parity violating models avoid metastability in the UFB direction? 

We will base our discussion on the Constrained MSSM (CMSSM) which contains the 
usual .R-parity invariant MSSM superpotential, 

W M ssm = h v QH 2 U c + h D QHiD c + h E LH 1 E c + nH x H 2 , (1) 

and a degenerate pattern of supersymmetry breaking with universal mass-squareds (171^), 
trilinear couplings (Aq) and gaugino masses (M1/2) at the GUT or Planck scale. The 
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notation is the same as in Ref . |3| . The extension to other patterns of symmetry breaking 
is straightforward. In addition we will consider only the two most important constraints 
on the CMSSM at low tan/? which are those due to Komatsu and Casas et al 0. 
(A similar analysis could be carried out on the other directions.) The corresponding flat 
directions we will refer to as directions 'K' and ! C respectively. 
Direction K is given by [0 

h° = -a 2 fi/h Dj 
d Lj = d Rj = ay./ h Dj 

v { = ay/l + a 2 y/h D] (2) 

for any pair of generations i, j, where a parameterises distance along the flat direction. 
(This does not give the 'fully optimised' condition || f|, but the difference will not be 
important here.) The potential along this direction is F and D-fiat, and depends only on 
the soft supersymmetry breaking terms; 

t r A 4 2/2/ 2 i 2 \ i 2 , 2 , 2 \ i o\ 

V = -tf~ a ( a ( m 2 + m Lj + m U + m d, + m Qj ) • ( 3 ) 

At large values of a, therefore, the potential appears to be unbounded unless m\ + m\. > 
0. Once radiative corrections are taken into account (and the very weak assumption 
that ml + rn 2 L . > at the GUT scale made), one instead finds that certain choices 
of soft supersymmetry breaking parameters cause the potential to develop a minimum 
radiatively at some high scale which depends on the weak scale parameters such as fi. In 
the CMSSM at the quasi-fixed point for example the resulting traditional UFB bound is 
m >M 1/2 H|. 

Direction C corresponds to, 

hi = -a 2 ^/h Ej 
£l 3 = e Rj = an/ h E] 

Vi = aVl + a 2 n/h E] (4) 

where i ^ j. The potential along this direction is 

V = ^a 2 (a 2 (m 2 2 + m 2 L .) + m 2 Li + m| +m 2 ). (5) 

In both of these cases we work in the basis in which the relevant Yukawa coupling is 
diagonal. The bound from direction C is comparable to that from K. 

In the next section we address the cosmological aspects of these flat directions, and 
argue that the answer to the first of our questions is no. That the UFB bounds are not 
necessary is already well known; the decay rate into a UFB direction at finite temperature 
is suppressed by a large, temperature-dependent, barrier in the effective potential 0. To 
show that it is not sufficient either, we need a more detailed understanding of the potential 
which we obtain through our analytic approximations. 

Replacing the traditional UFB bound with a necessary condition would require de- 
tailed knowledge of early cosmology which is, unfortunately, lacking. We shall instead 



3 



propose a truly sufficient condition which arises from requiring that there is only one 
minimum in the zero-temperature effective potential and that it is at the origin. By 'suf- 
ficient', we mean that we do not have to appeal to any cosmological model (even the so 
called Standard one) to end up in the correct minimum. 

In the third section we develop an analytic treatment of both the traditional and 
sufficient bounds. We explain why they are numerically very close (and can in fact only 
be resolved with two-loop or better accuracy). 

In the fourth section we address the second of our questions. It arises because, as is the 
case for nearly all F and D-flat directions, the K and C directions correspond to (analytic) 
gauge-invariant polynomials which are absent from the superpotential 0. Stated in this 
way, it is clear that it is discrete symmetries that lead to UFB problems/metastability. 
Of course some discrete symmetry is required to prevent baryons decaying. However, in 
models which violate i?-parity, the most dangerous directions can be lifted, although it is 
not clear a priori whether it is possible to add all the required couplings without violating 
experimental bounds. It turns out that the necessary couplings can be added and destroy 
the dangerous minima if the i?-parity violation is just below experimental bounds. 

2 Vacuum decay into and out of flat directions 

Before getting to grips with detailed calculations of the bounds, we will address the first 
question above which concerned the relevance of the traditional UFB bound. To do this 
we need to find the decay rates into and out of the dangerous K and C directions. In 
this section we estimate the false vacuum decay rate into the radiatively induced minima 
at finite temperature. We also consider what happens when the new minimum which 
develops radiatively along the flat direction is not the global minimum. In this case 
we find that the decay rate back to the global minimum at the origin is also greatly 
suppressed. Hence we will conclude that whether the new minimum is global or local 
(the criterion of the traditional UFB bound) is irrelevant. The only relevant bound we 
can impose (apart from the need to avoid quantum tunneling into the charge and colour 
breaking minimum which only excludes very low values of m [5] ) is a sufficient one - that 
there be no minima at all along the UFB directions except the physical one. 

We first borrow some results from later sections to get some idea of the shape of the 
potential along the F and D flat direction. Figure 1 shows 

V = Vhl/fMl /2 (6) 

in the K direction (with j = 3) in the CMSSM for the three values of mo indicated. It is 
plotted in the approximation that we can neglect the bottom and tau Yukawa couplings 
and two-loop corrections in the running of parameters. (A full two-loop numerical running 
produces something similar.) We have taken fi = 500 GeV, mi,(m t ) = 2.5 GeV and 
ot3(mt) = 0.108. The potential is taken at the quasi-fixed point so we don't have to 
worry about the value of Aq but is otherwise typical. As can be seen, a new minimum 
forms radiatively when mo/Mi/ 2 ^1.07 and becomes global when mo/Mi/ 2 £ 1.05. The 
latter value corresponds to the traditional UFB bound (modulo the approximations we are 
making here). The minimum in the flat direction (which we shall denote by 0) forms when 
a = a m i n = 8.5 so that <p m in = few x 10 6 GeV. One fact which will prove extremely useful 
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when deriving the bounds analytically later is that it is the first term in the potential 
(i.e. m\ + m 2 L ) which dictates the form and position of the radiatively induced minimum. 
The depth of the minimum is therefore ~ n^w^min = a^^m^/h 2 . 

Now consider the decay rate into the new minimum at finite temperature. In Ref. U it 
was argued that the traditional UFB bounds are not necessary because the decay rate into 
the UFB directions is generally far too small. The suppression occurs due to the presence 
of a large, temperature (T) dependent barrier between the false vacuum (where we are 
living if the UFB bound is violated) and the true charge and/or colour breaking one. This 
barrier is generically present due to the flat direction, 0, giving masses to fields to which 
it couples. For 0^,T the contribution to the temperature dependent effective potential 
from a field, i, is suppressed as e~ mi ^l T (where m; is the mass it derives from coupling 
to the flat direction) and the potential is the same as that at zero temperature with a 
curvature of order rriyy coming from the soft supersymmetry breaking terms (c.f. Eqn.(0)). 
For ~ the contribution to the effective potential is 0(— T 4 ). The escape point of the 
critical bubble is therefore of order ~ T 2 /mw Since the barrier height is of order 
T 4 the decay rate is suppressed by a factor e~ S3 ^ T where S3 ~ T 4 /m^, which is huge 
unless T < mw- At zero temperature the quantum tunneling rate is only large enough for 
vacuum decay to have taken place within the lifetime of the universe for very low values 
of m f|. 

The estimates of Ref.|| apply equally to the dangerous K and C directions we are 
considering here. Actually we can make a much more accurate estimate of the action 
for the finite temperature decay when < T < 4> m in as follows. The one-loop finite 
temperature contribution to the effective potential is 



V T (ct>) 



rp4 

2^ 
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1 =F exp (- v /g 2 + mf(0)/T 2 



(7) 



where the upper sign is for bosons, the lower one for fermions, and the rij are the cor- 
responding degrees of freedom. The masses are field dependent parameters which 
generally have a wide range of values and increase roughly linearly with distance along 
the flat direction. The potential to calculate the decay rate is normalised to zero at the 
origin, 

AVH0) = V T {<f>) - V T (0). (8) 
The integral has the well known limits when mj C T and rrii 3> T; 



AV T (cf>) 



^ E [n B Am 2 B + ?fAm 2 F \ ; m,(0) < T 
vr 2 ^T 4 /90 ; m,(0)>T 



(9) 



where Am 2 = m 2 (0) — m 2 (0) and 



g* = n B + (7/8)n F 



(10) 



counts the difference between the degrees of freedom of the heavy particles at that partic- 
ular value of (i.e. those for which m« ^> T) and at = 0. Particles which are strongly 
coupled to the flat direction cause a steep rise in the effective potential near the origin. 



The contribution near = can also be expressed as V(4>) 



rn; 



: ,)0 2 /2 where 
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and 



is the (negative) zero temperature mass-squared of the F and D-flat direction 



m 7 



24 



E 



»bA b + — A F 



in which Xb,f stands for the coupling of the field to the flat direction (~ g 3 , gi or h t in this 
case). (We take the mass-squared to be constant (and always negative) for this estimate.) 
The finite-temperature contribution to the effective potential of particles which couple to 
the flat direction can therefore be approximated by 



V&) = -<P 2 T 2 Q{kT - k'<P) + x 9 



where we have defined 



K 
. '2 
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K - ml)/T 2 
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2 e(«v-KT), 



(12) 



Trg*/45. 



(13) 



The first term is the T limit, the second term is the <p^> T limit, and the last term 
is the zero temperature contribution which is still present when 3> T. (The second T 2 
term is just to match at k'(/) = kT.) For the moment we are assuming that all the particles 
couple to equally and one should bear in mind that in the MSSM there are of course 
many different couplings. To get the decay rate we find the 3-dimensional euclidean action 
for a (spherically symmetric) critical bubble, 0(p) where p is the radial coordinate. The 
bubble is a solution of the equation M 
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2<90 
pop 



dV 



with the boundary conditions 



(14) 



oo 



dp 



0. 



By matching and d(j)/dp at k'0 = nT, we find the solution 

{kT p sinpm^ ~ 
k' p sin pm^ " ' 

where the radius of the bubble is given by 



p= — f 7r — tan 1 (m ( h/KT)) « ir/mj, 
m v 7 



and the escape point by 



0(0) 



TTK 2 T 2 



pm^nT 
k' sinlpm^) k' m$ 



(15) 



(16) 



(17) 
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where the approximations are for kT 3> m^. The three dimensional euclidean action is 
then given by 

S 3 /T = 4./^(lg) 2 + ^) 



2ttV2 ^3 
3 m 3 
2tt 6 ^ T 3 



135 m 



3 ' 



(19) 



This action is independent of the precise couplings to the flat direction because we assumed 
that the temperature is high enough for it to be dominated by the very large escape point. 
Hence the action is the same when there are many different couplings to the flat direction, 
Aj, provided that they are all either strong enough for particle i to freeze out well before 
the escape point or weak enough for it to not contribute appreciably to the effective 
potential. These two conditions can be written; 

either Aj ^> m^/T 

or Aj < rn<p/<f) min . (20) 

The criterion for the vacuum to have decayed within the lifetime of the universe is 
(see for Ref. f| for example) 

S 3 /T< 200. (21) 

For T 3> the action is indeed much more than 200T and becomes larger for higher 
temperatures, from which we conclude that the false vacuum is unlikely to have decayed 
within the lifetime of the universe. At T ~ we would need to take into account all 
the field directions in order to estimate the rate as well as the depth of the electroweak 
symmetry breaking minimum which we were able to neglect above. However from the 
above and Ref. f| it seems unlikely that the false vacuum would decay within the lifetime 
of the universe even in this case. 

Now consider decay from the new minimum to the origin. This would be relevant if 
either the new minimum were local {e.g. if 1.05Mi/2 < mo < 1.07Mi/2 in figure 1) or if 
it were global at zero temperature but lifted by finite temperature corrections (of order 
T 4 where T <C <fi m in) ■ The barrier height (in figure 1 for example) is governed by the first 
term in the zero temperature potential, and is therefore of order af nin p 2 my V /hl, and its 
width is of order min = a^ in /i//i{,. The critical bubble is determined by Eq.(|TJ) in which 
the RHS is of order (a^^p/h^m^. We can obtain an estimate of the action by scaling 
out the dependence on parameters using the dimensionless variables (p = 4>hb/ (a^^p) 
and p = my/p. At finite temperature a critical bubble has a three dimensional euclidean 
action given by 



a 4 u 2 



S 3 = 47r-f^/ (22) 



where / is the integral, 




LL minP l,l W 
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which we expect to be of order 1. The decay rate is therefore only significant for large 
temperatures. Decay out of the false vacuum (for the potential of figure 1) within the 
lifetime of the universe would require, 

4 2 

200T > 4tt^P^ ~ 5 x 10 11 GeV. (24) 
h b m w 

Of course our approximation breaks down long before this temperature is reached since, 
when T ^> <f) m i n ~ 10 6 GeV, the finite temperature contributions to the effective potential 
(proportional to <p 2 T 2 for T ^> <p) will in any case destabilise the local minimum. The 
quantum tunneling rate is also suppressed; for significant tunneling within the lifetime 
of the universe the 0(4) symmetric 'bounce' solution || should have an action S4 <,400. 
We estimate 2 

S 4 = 27r 2 ^f^-/' > 400 (25) 
h b m w 



where the integral 



/o I 2 \dp a^i^mfy 



which is again expected to be ~ 1. 

With these two pieces of information to hand, it is clear that the traditional UFB 
bound is irrelevant; it can never be either necessary or sufficient. It is not necessary 
because if the universe were trapped in a metastable minimum at the origin it would not 
have decayed within its lifetime. It is not sufficient because had the universe been trapped 
in a radiatively induced local minimum it would also not have been able to decay. 

Finally, we can address one of the suggestions in Ref.|| for driving the universe to a 
metastable vacuum at the origin. Some sort of dynamical mechanism is probably required 
if such a scenario is to avoid fine tuning. This mechanism could simply be a sufficiently 
high temperature, Tr, in the early universe. Each field, i, then contributes positive curva- 
ture to the potential for Aj0 < Tr. (Other proposed mechanisms depend on the additional 
soft supersymmetry breaking terms which can be generated during inflation |10|.) How- 



ever the destabilisation of the radiatively induced minimum can only be guaranteed if T R 
is greater than the zero-temperature position of the radiatively induced minimum, ^) min . 
For the examples we are considering, 

T R > few x 10 6 GeV. (27) 

However when, as we are assuming, supersymmetry is communicated gravitationally there 
is a bound on the reheat temperature coming from nucleosynthesis; if Tr is too high, 
nucleosynthesis is disrupted by the creation and subsequent late decay of gravitinos 
This translates into the bound, 

T R < 10 9 GeV. (28) 

Thus, when mo is close to the traditional UFB bound, the reheat temperature has to fall 
inside the window which is left. Moreover, this window closes when the traditional UFB 
bound is more strongly violated. Figure 2, which we will derive shortly, shows the position 
of extrema along direction K at the quasi-fixed point plotted against m^/M 2 ^ (with the 
same choice of parameters as in figure 1). To the right of the cross are only minima, to the 
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left maxima. The cross corresponds to the point of inflexion in the mo = 1.14Mi/2 curve 
of figure 1. As m\ decreases below the traditional UFB bound, the minimum deepens 
and (f) m in grows rapidly. For low enough uiq, Tr has to violate the bound in Eq.(p8|) if 
it is to destabilise the radiatively induced minimum. (For the quasi-fixed point example 
in figure 2, the window disappears when m < 0.7M^ 2 .) If the reheat temperature is far 
below (pmin, then we have to ensure that the contribution to the effective potential from 
the more weakly coupled fields (proportional to \ 2 cj) 2 T 2 with A, <C 1) is enough to drive 
(p to the origin, a complicated and model-dependent task. 

To summarize, without a detailed knowledge of early cosmology, there is no well- 
defined criterion for excluding regions of parameter space based on the existence of deep 
minima along F and D-flat directions, apart from at very small values of where 
quantum tunneling can take place within the lifetime of the universe. In the following 
sections therefore, as well as deriving the traditional UFB bound analytically, we shall 
focus on a sufficient condition which is better defined, namely that there are no minima 
other than the one at the origin. This condition ensures that the origin is reached from 
any initial even at zero temperature. The areas of parameter space which satisfy this 
condition are therefore allowed regardless of cosmology, and figure 1 suggests that, at 
least for the UFB directions we are considering here, they may not be very different 
from those allowed by the traditional UFB bound. At the quasi-fixed point for example, 
figure 1 gives the traditional bound to be m •> 1.05Mi/ 2 and the sufficient condition to 
be m •> 1.07Mi/2. Thus, in the one- loop approximation we are making (and possibly to 
two- loop order as well), the traditional UFB bound is indistinguishable from the sufficient 
condition. 

3 Analytic derivation of 'UFB' bounds 

In this section we develop an analytic treatment of both the traditional bound and suffi- 
cient condition coming from the K and C directions in the CMSSM at low tan (3. These 
are the two most restrictive constraints and were analysed numerically in Ref. ||. We 
first summarize the one-loop RGE solutions, and then use them to obtain the bounds. 

We shall obtain the bounds two ways in an order which may at first seem a little odd; 
we first present an accurate and straightforward but semi-graphical method, and then we 
present an entirely analytic method that requires some approximations. The reason for 
this ordering (it is after all more usual to first derive something approximately and then 
do it accurately) is that the second method is far more powerful, allowing us to produce 
bounds which are valid for all low tan (3. It also gives us more insight into how the bounds 
depend on the unknown parameters (such as //) and why they are numerically so close to 
each other. In fact we will conclude that (as anticipated in the previous section) two-loop 
accuracy is required to resolve the traditional bound and the sufficient condition. 

To be able to do an analytic treatment, we will need to express the mass-squared 
parameters appearing in the potential as a function of scale {i.e. not as a numerical 
integral). This is possible if we work to one- loop and in the approximation that we can 
neglect couplings other than a^j «3 and h t in the running. We first give solutions for all 
the parameters in terms of the three quantities which are strongly attracted to quasi-fixed 
points in the minimal supersymmetric standard model (MSSM) (although other terms 
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formally have quasi- fixed points too). The three parameters with quasi- fixed points are 

R = h 2 Jgl 
A t = A Us3 

3M 2 = m 2 U33 +m 2 Q33 +ml (29) 

When h t is high at the GUT or Planck scale (higher than 1 say), these three parameters 
are completely determined at the weak scale irrespective of the pattern of supersymmetry 
breaking at the GUT scale; 

E? FP = 0.87 

A QFP = _ lmMl/2 

(M 2 ) QFP = lMM 2 /2 . (30) 



Their values govern the running of the MSSM at low tan/? [I2|. (As pointed out in 



Refs. fl3|, |6| and references therein, various constrained versions of the MSSM lead to 



many more quasi-fixed points in addition to these three ubiquitous ones, and can lead to 
enhanced predictivity and reduced FCNCs at low tan/3.) 

For completeness we present all the solutions for the flavour diagonal parameters here 
and for the moment leave in the ai dependence; defining 



ai = r 



di\o 

g = ^3 43/9 « 2 - 3 «r 5/99 (si) 



where the 0-subscript indicates values at the GUT scale, these are 

A Utl - U, = M 1/2 (-?«« + §4" + + (Au„ - \A 



t \0 



I . / 40- 5„ (1) 29„ (1) \ , . 1 



Ab - -A t = M 1/2 (--5^ + -5^ + -5^)+(A b --At)\o 

A Dii = Mv^-^ + S^ + ^S^+^aJIo 

A Eaa = M 1/2 (38? + ^S{^ + (A Eaa )\ 

B-\At = M^^ + ^ + ^+^-^lo 

< 33 ~M 2 = M^(^ 2) + ^-^) + (< 3 -M 2 )| 



m?>„. - -M 



2 

-^M 2 = M 2 /2 (-^ 2) + ^! 2) ) + (m 2 -^M 2 )| 
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m 2 




m u 


= M lnhf- 

\ y 


m o 




m D 




n 





^ 2, J + «,)lo 
2 

99 



*i 2) ) + «,jlo 



(32) 



where i — 1,2 and a = 1, 2, 3. 

The solutions of the remaining three parameters R, A t and M 2 , can be expressed as 
functions of 

r=l/a s = l + ^lo gx ^-, (33) 

so that 

3 ~ 5 

«2 = -j ; «i = — R — — • 34 

4 — r 16 — llr 

Taking a^{m t ) = 0.108 means that 0.37 < r < 1 with r = 1 corresponding to the GUT 
scale. We shall use r' 1 and d 3 intechangeably. If we define 

tt/ \ ~ 16/9- -3 ~ -13/99 

ll(r) = a 3 a 2 a x 

J = — I— f\{r')&r' (35) 



r n(r) 

then, in terms of the GUT scale values (R , A , m 2 ,) plus the functions 

R QFP = 1/2 J 
R 



P 



R QFP 

1 — r 

t = — p-l (36) 



we find 



1 1 1 

+ 



R RQ pp R Ilr 

A = (1 - p)A + Mi /2 (p£ - (1 - r) (—& 3 + «2 + ^r&i 
M 2 = (l-p^-^+^Mf/^l-rh (37) 



where 



K = (l-p)(eM V2 -Ao) 2 -e 2 M 1 2 / 2 + M 1 2 / 2 (l-r)(e + l)(y<5 3 + a 2 + ^<5 1 ) 

16 13 

7 = a 3 (l + a 3 (l-r)/2) + a 2 (l-a 2 (l-r)/6) + — a 1 (l-lla 1 (l-r)/10)(38) 

y 4o 



ll 



There is only one integral to do, J, which we must now approximate by dropping the at 
terms 

16479 r" 7//9 

Jrn « - — — + — — (14080 + 36960r - 1680r 2 + 77r 3 ). (39) 
1540 4620 

This approximates the full integral to about 5% in the region of interest (i.e. where the 
minimum develops). 

To make contact with the usual parameterisation of the CMSSM we can express the 
GUT scale parameter Rq as a function of of the weak scale parameter tan f3 for which we 
need to insert the running top quark mass, m t (m t ) = 167 GeV, and use 

R = A m 2-*« ( 4 °) 
A7ra 3 v z sin p 

where v = 174.1 GeV. Inserting the solution for R gives 

p-i = n^^l^fsin^-sin 2 ^) 

= 6.34 (sin 2 p- sin 2 (3 QFP ) . (41) 



In this approximation (tan/3)^ FP = 1.89. However the value we infer for Rq is clearly 
highly dependent on its exact value and hence is sensitive to two-loop corrections. In 
this sense Rq is a more natural parameter to use than tan/3; the latter is equivalent to a 
very peculiar weighting of the GUT scale parameter, so that mid-range values of tan f3 all 
correspond to roughly the same minimum value of Rq. As a nod to convention, however, 
we will present some of our results in terms of tan /3. 



Now let us use these expressions to examine the potential of Eq. (0) in more detail for 
the case of the CMSSM in the j = 3 direction. We first use a semi-graphical method to 
extract the traditional bound and sufficient condition. 

The mass-squared parameters run with scale, and the appropriate scale, A, to evaluate 
them at is where the finite one-loop contributions (which we are neglecting) to the effective 
potential are small, A ~ h t <f) (assuming that the minimum occurs for a ^> 1 which we 
shall verify in a moment). This should correspond roughly to the largest mass if we were 
to diagonalise the mass matrices and evaluate the one-loop contribution completely, a 
can then be expressed as a function of r, 



a 2 (r) = b(r) exp(c(r - 1)) (42) 

where 

b(r) 

c = 47r/6a 3 |o = 52.4 (43) 



M GUT h b _ 2.85 x 10 12 / Mqut \ ( m b (r) \ /200GeV\ 
h t fi ~ h t (r) /l.l V2 x 10 16 GeV/ V2.5GeVyl V /x(r) j 



where the central values are those we shall adopt for the scale A = m t or r = r t . We can 
solve for b(r) in the same manner as above, and find 

1/ \ ,i / R \ 2 ^ ~ 14/27 ~ -1/2- -7/594 

K r ) = b \o (j^) a 3 ' a 2 a x 
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{t) \R(r t )) \a 3 (r t )a 2 (r) ) 

Denoting the position of the radiative minimum by <f) m in, 
saturated by 

dV 



1/2 

(44) 

the traditional bound is 
(45) 



d 2 V 



(46) 



and the sufficient condition is given by 

dV_ 

and so we can work with 

V = Vh 2 h fo 2 M 2 1/2 (47) 

which depends only on b(r). (This means that for the sufficient condition we are neglecting 
the running of h t which is relatively small in this region.) The potential becomes 



V(r, Too) 



a 2 (r) 



Ml 



a\r)C-M 2 + -mi - M* /2 (^f + ^ + — 



f( 2 ) 



/2 

-M 2 



-5^ 



29 
594 



(2)n 



(4S 



Figure 1 shows V for m = 1.03Mi/ 2 , 1.05Mi/ 2 , and 1.07Mi/ 2 , for the central values 
of parameters we used in Eqn.(|43|) and for the quasi- fixed point where Rq 1 = and 
where there is no dependence on A in the mass-squareds. The minima of V are at 
(fimin = few x 10 6 GeV. At the quasi-fixed point the traditional bound is saturated with a 
minimum corresponding to a m i n = a(0.5427) = 8.5, justifying our assumption of a 2 ^> 1 
in direction K. 

As can be seen, for these values of m the minimum is being rapidly lifted, with 
m = 1.05Mi/ 2 saturating the traditional UFB bound. From this figure we have already 
anticipated one of our main conclusions; radiatively induced minima along UFB directions 
are lifted very quickly by changes in the parameters (mo in this case). This happens 
because the negative a(r) 4 term in the potential balances the positive a(r) 2 one when the 
traditional UFB bound is saturated, so that the height of the minimum is varying rapidly. 

Ideally one would like to be able to locate the extrema for a given mo, Ro and Aq, 
see for what value of mo there exists a minimum away from the origin, and hence derive 
the traditional UFB bound. To make the problem tractable analytically we can use 
the fact that there is a one-to-two correspondence between m and the location of the 
extrema, r = r min . Consider again the K direction, with a minimum at 4> m i 
corresponding to 



a minf J ' fob 



6a(M t 



GUT) 



47T 



log 



(49) 



and denote the value of m corresponding to an extremum at r = r min by m (r m 
Solving for ml(r min ) is relatively simple since V(r,m ) depends only linearly upon it; 



"'o(''™i) 



Ml 



/2 



d(a 2 + 



(a 2 + \)R{r r 



r^ n (3M 2 + A; 



d{a 2 + -)M 2 
6 
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+ — a\ + 2(a 2 + l)a\ + (la 2 + U ] :! 



,5 15, 



2. 



where d(r) encompasses the running of a 

1 da 2 dlogb ^ 

a 2 dr dr 

1 1 4 

= -« 2 + -oc\ - —R + c. (51) 

2 6 3r 

Eq.(|50D is implicit in that, away from the fixed point, M 2 itself depends upon m\. How- 
ever, since the dependence is (by dimensions) also linear the full (but unwieldy) expression 
can easily be obtained from the above. 

For a given value of Aq and Rq (or equivalently Aq and tan/3) we can now find the 
sufficient condition from Eq. (p0|) . First consider the quasi-fixed point. Figure 2 shows 
™o( r m»)M 2 /2 plotted against r min . Values of r min to the left of the cross are maxima, to 
the right, minima of the potential for a given value of m . This function has a maximum 
when m\ = 1.14M 2 , 2 so that the only minimum when mo > 1.07Mi/2 is at the origin, 
as we saw in figure 1. Close to the quasi- fixed point the sufficient condition can be 
approximated by expanding in 1/Rq and by noting that, from the analytic solutions, they 
can be expressed as 

4 > ^ + ( c ' + ^ +c 4) / ^ +0(lK) ' l52) 

By matching at the quasi-fixed point and at Rq = 5 and A = (—1, 0, 1), we find that 
c . .3= 1.14, -1.69, -0.63, 0.30 respectively. These numbers approximate the bound (in the 
one-loop approximation) to better than 5% for i? ~4 (corresponding to h t {MavT) > 1-4) 
and \Aq\ <Mx/2 and of course much better for larger R . We can also find an 'absolute' 
condition by completing the square in Eq. (|52|) ; values of m which satisfy the sufficient 
condition in direction K must also satisfy 

2 

— f - > c + (cx - c 2 2 /Ac 3 )/R = 1.1 - 2.0/^0- (53) 

M l/2 

If 777.Q is below this bound a minimum is radiatively generated whatever the value of A . 

We can now find the traditional UFB bound graphically by plotting the height of the 
extrema, V (r min , m (r min )) , as shown in figure 3. At the quasi-fixed point, the minimum 
becomes negative at r min = 0.5427, giving the bound m ^,m (0.5427) = 1.05Mx/ 2 - Again 
by matching at the quasi-fixed point and at i?o = 5 and Aq = (—1, 0, 1), we find the 
coefficients Co.. 3= 1.10, -1.67, -0.63, 0.30. Thus the sufficent condition is as expected 
virtually indistinguishable from the traditional bound. 

Applying the same treatment to direction C at low tan (3 (which was studied numeri- 
cally in Ref.0), we find co..3= 1.15, -1.54, -0.57, 0.26 for the traditional UFB bound and 
for the sufficient condition we find Co.. 3= 1.17, -1.56, -0.57, 0.26. (For this direction the 
definition of a 2 the same, but with b(r t ) = 1.9 x 10 12 for m T (r t ) = 1.7 GeV.) These bounds 
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are in agreement with Ref . given our neglect of two-loop contributions to the running of 
the gauge and Yukawa couplings. Ref.[|J finds m £Ml/ 2 at the quasi-fixed point (which 
should be compared with our result above of m > 1.05Mi/ 2 )- 

Finally let us consider the accuracy of the approximations we made. We saw that 
radiatively induced minima along UFB directions are lifted very quickly by changes in 
m . This means that the bounds are less sensitive to our lack of knowledge about the 
parameter b and in fact depend on it only logarithmically as we shall see shortly. (In 
fact the error involved is of the same order as that we associate with neglecting the finite 
one-loop corrections.) We can also see that the bound increases with b as we would expect 
since increasing b increases the negative contribution to the potential from rnj. This is 
in accord with Ref.Q where the effect of increasing the 'unification' scale to the Planck 
scale was examined. We also see that the bound increases with mj (which is why we chose 
j = 3) and decreases with \i. The variation of the bound with A is quadratic becoming 
stronger at large tan/5 with a minimum close to A = at positive A . 

Although the above method is accurate for a given choice of A and tan (3 (given the 
other approximations we have made), it is cumbersome and in particular is not well suited 
to finding bounds at smaller values of Rq. In addition we would like to understand why 
the sufficient condition is numerically so close to the traditional bound. Because of this 
we now turn to a completely analytic derivation of the same bounds. 

For this we can use the fact that the minimum appears at a scale where the first 
term in the potential is driven negative by m|. So, considering direction K, if we define 
A(r) = (m| + fn 2 L z ) / Mlj 2 and B(r) = (m 2 L . + m 2 d . + mQ.)/My 2 , the minimum of the 
potential is at a scale close to the point which we shall call r p where A becomes negative, 
A(r p ) = 0. We define two additional scales; = the position of the minimum when 
the traditional UFB bound is saturated (i.e. when the minimum is degenerate with the 
origin), and r c = the position of the point of inflexion when the sufficient condition is 
satisfied. We first determine these scales in terms of r p as follows. It is convenient to work 
with 

V = M2 % 2 V = (e^A + e^B/b) , (54) 

m 1/2 1V1 GUT 

(Note that instead of Mqut we could have used any scale, Ax, as long as c = 4tc /Qa 3 \\ x .) 
Since the scales r p will appear algebraically in the determination of the bounds we can 
sacrifice a little accuracy in their determination, and approximate the potential close to 
r p by developing A(r); 

V(r) « (r - r p )A{r p )e 2c{r - 1] + e c(r " 1) 'B(r p ) /b(r p ). (55) 

Because at this point A is running much faster than B, b and h t , we are going to neglect 
the running of B, b and h t close to r p . For the bounds we will need 

V(r) w (l + 2c(r -r p ))Ae 2c(r - 1) +e c{r - 1) cB/b 

V{r) ss Ac{l + c{r -r p ))Ae 2c{r - 1) +e c{r - 1) c 2 B/b. (56) 
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For the traditional UFB bound, (V(rd) = V(rd) = 0), we find 

, 1, fecB\ 
r p = H — m 



c \Ab) 

r d = r p -l/c. (57) 



Likewise, for the sufficient condition, {V{r c ) = V(r c ) = 0), we find 

1 1 (e^ 2 cB\ 
r v = H — in 



c \ 2Ab 



ib J 



r c = r p -3/2c. (58) 

Since c is large (c.f. Eq . (f43"D ) the minima are very close to r p (as was the case for the 
example shown in figure 1). The value of r p in each case is determined using the expression 
for B above and for A which follows straightforwardly from the renormalisation group 



equations [ITJ 



A = -(3M 2 + AT) - 2a% - -a\. (59) 
r ' 5 



A comparison of the expressions for B and A together with A(r p ) = gives 

A/B » R/r (60) 

to a good approximation independently of mo/Mi/ 2 . This allows for a sufficiently accurate 
determination of r p as the dependence on the parameters above is only logarithmic. For 
relatively large Rq £ 1 it is well approximated by using the quasi-fixed values of parameters 
in the logarithm 

nmm h K 7 2 , A/r 2 \ , / 0.544 traditional UFB , , 

r p = 0.019 log(1.67 + m /M 1/2 ) + j Q ^ (61) 

The values of r„ are those which must satisfy 



and solving this gives the bounds on we are looking for. We find the approximate 
expression for the bound 



4p p - p 2 - 1.67 + p p (l - p p )a 



Mf /2 4 - 3p p 

where 

1 i?« FP 
1 + 



(63) 



Pp .Rollr 

% Ml/2 



3i?n 



« ~ tt^-2tt^- (64) 
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We again see from Eq.flB3|) that the bounds are always quadratic in A (except for the 
very weak A dependence in the determination of r p ). We also see that the bounds are 
much more restrictive for negative Aq and have a minimum at A rs M1/2. 

As an example of the usefulness of this method we show in figure 4 the 'absolute' 
bound on m in this direction for all tan/3 < 15 where our RGE solutions are accurate 
to one-loop order. For values of m below this bound there is a radiatively induced 
minimum whatever the value of A . Replacing p p in Eq . (|63"D with tan/3 gives a nice 
approximation to this diagram. (Since the diagram goes well away from the quasi-fixed 
point we found r p by solving the first equations of Eqs.( |57| , |58D fully rather than using the 
large Rq approximation above.) The plateau in figure 4 is a result of the mid-range values 
of tan/3 all corresponding to roughly the same value of Rq. 

To compare this approximation with the previous method, we find the coefficients for 
the expansion near the quasi-fixed point of c ..3= 1.13, -1.67, -0.63, 0.30 and c ..3= 1.09, - 
1.65, -0.63, 0.30 for the sufficient condition and traditional UFB bounds respectively. 

We can now see why the traditional UFB bound and the sufficient condition are so 
close. The separation between the scales and r c at which the minima form can only be of 
order 0(1/ c). When we solve for A(r p ) = the values of r p differ by a factor e 1//2 /2 = 0.82 
appearing in the logarithm. Thus we would have to know b and in particular /i to better 
than 18% (i.e. to two-loops) before we could distinguish the two bounds. On the other 
hand the parameter r p depends only logarithmically on b but appears algebraically when 
we determine the bounds - hence the bounds on mo can only depend logarithmically on 
b. 



4 Lifting flat directions with imparity violation. 

In this section we will see how much the flat directions can be lifted if we choose 'baryon 
parity' rather than i?-parity to prevent proton decay JTH . 

As we stated in the introduction, violating i?-parity can lessen UFB problems because 
there are less gauge invariants. The role of gauge invariants in flat directions was discussed 
in Ref . 0] . To clarify the role of the i?-parity violating couplings let us briefly summarize 
that work. If we represent the scalar component of a generic superfield, k, as z^, then the 
potential contains D 2 and F 2 terms where 

D A = z k T k A lZl (65) 

and T A are the matrices of the generators of the group. In addition, gauge invariance of 
any analytic function I(zk) requires that 

hT A zi = (66) 

where the subscript on the I implies differentiation by Z)~. Thus, there exists a D-flat 
direction corresponding to every 

h\(z) = G{z k ), (67) 
where C is a complex constant. The inverse, i.e. that every flat direction corresponds to 



an invariant satisfying Eq. (67), is a less trivial property. If the invariants in question do 
not appear in the superpotential, then the F 2 can often be made to cancel. 
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If we didn't already know the UFB directions in the MSSM, then we could have 
used Eq.(|67D to do it methodically (which is extremely useful for more complicated flat 
directions). For the K directions, the relevant invariants are 

W = H x (h D] QjDj + fxH 2 ) 
I = L i {Q j D j + pL'H 2 ) (68) 

where W is a term in the superpotential while / is not allowed if one assumes -R-parity. 
The K direction in Eq. (12) is found as the solution of Eq. fl6T|) for / with the coefficient // 
being chosen to cancel F Hl , ji 1 = —a^fi/ho- Similarly, the C direction can be constructed 
from the -R-parity violating invariants 

LiLjEj ; H 2 L { (69) 

with i 7^ j. 

Notice that the absence of a particular gauge invariant operator from the superpoten- 
tial is not enough per se to ensure F-flatness unless no pair of superfields in the gauge 
invariant operator appear together in any term of the superpotential. For example Qj and 
Dj both appear together in the Yukawa couplings of the MSSM and this is why we had 
to construct our flat direction by the conjunction of two invariants in order to cancel F^. 
An example of a single gauge invariant which corresponds to an F and .D-flat direction is 
the other .R-parity violating operator, 

eUDD. (70) 

This flat direction (although it has been discussed at some length in the literature) is 
usually safe in models of gravitationally communicated supersymmetry breaking. 

Baryon parity allows all of the invariant operators above corresponding to the K and 
C directions, and so can potentially lift all the flat directions we have been discussing. 
Since the K and C directions with j '• = 3 result in comparable bounds on parameter space, 
we will assume that all the directions should be lifted in order to significantly change 
the conditions on parameters seen for example in figure 4. There are 5 flat directions in 
total, corresponding to i = 1, 2, 3 in the K directions and % = 1, 2 in the C directions. In 
this section we will attempt to lift all five, whilst satisfying experimental constraints on 
.R-parity violating operators. 

We begin by adding to the superpotential the additional operators 

Wb = \jkLiLjE k + X' ijk LiQjD k . (71) 

We do not consider the i?-parity violating bilinear operators, /ijLjif 2 , since the dimension- 
ed coefficient /ij must be much less than mw and so their contribution to the potential 
at large field values will be very small. Along the K and C directions the contributions 
to the potential from the F 2 terms are 

v F - terms = E^(i^33r+(i+a 2 )iAL 3 r+(i+a 2 )iA: 3m r) 

m "'6 

4 4 

V F -ter m s = £ 2^ (| A m33 1 2 + (1 + a 2 ) \ X im3 1 2 + (1 + a 2 ) | X l3m | 2 ) (72) 

m t 
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respectively. In addition one can expect contributions from new .R-parity violating trilin- 
ear terms. Since the sign of the field (i.e. the sign of a) has not yet been determined, the 
potential will be minimised when they are negative; 



VA-terms ~ ~ I^A'A^/i 



\<i\"yl + a 



3 |a| 3 Vl + a* 

VA-terms ~ — \A\XiS 3 fI \ — (73) 

respectively. The F-terms must overcome both these terms and the radiatively induced 
minimum at m j n ~ a£, in ////i& if they are successfully to lift the flat direction. 

We now discuss how big the new couplings have to be, generically, to lift the flat 
directions. Let us begin by making a rough estimate in which we say that the new 
F-parity violating F-terms have to 'fill' the radiatively induced minimum; i.e. we are 
neglecting the fact that the new terms change the position of the minimum. (We shall 
see shortly that this estimate is surprisingly good.) Since the depth of the minimum is 
~ (imin^ 2 m w /titi the i?-parity violating couplings (generically denoted by A and A') must 
satisfy 

A>^ - 0.004^ 

y>W^ 0004 ^ (74) 

Q>minft ft 

The couplings A and A' are understood to be evaluated at the scale of the minimum. If 
some of the flat directions are lifted by the a 4 terms in Eq. fl72|) then the A and A' must be 
on the edge of experimental detection unless /i is very large [|l|| . We therefore concentrate 
on the a 6 terms and the condition above. Assuming that, as in canonical gravitationally 
communicated supersymmetry breaking, A\ ~ Ay ~ mw and that all of the -R-parity 
violating couplings the same order of magnitude, the trilinear terms are dominated by 
the F-terms if Eq. (|74j) is satisified and a min 3> 1. 

We can refine our estimate substantially by adopting the approach of the last section. 
If .R-parity violating terms are to make flat directions safe then we should, to be consistent, 
insist that there are no minima at all except the Standard-Model-like minimum for any 
value of parameters. Defining 

(75) 

(and similar for A') the form of the potential is now 

V = Vh 2 T /n 2 Mt /2 = a 4 A + a 2 B + a 6 A 2 , (76) 

where a(r), A(r) and B(r) are as before. We now find the sufficient condition - i.e. a 
condition which ensures that the F-parity violation has destroyed all minima along the 
flat direction. The sufficient condition is then 

=0 (77) 

da 2 da 
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where A and B are to be evaluated at the scale r c which we have yet to find. The 
determination of r c is simplified once we realise that all three terms are essential in the 
lifting of the minimum. The a 6 A term simply pushes the (extremely deep) minimum in 
the potential to low energy scales until it runs up against something positive - in this 
case the a 2 B term. Thus a suitable criterion for the determination of r c is simply that 
a 2 A(r c ) + B(r c ) = 0. Once we have r c we solve Eq.flTT|) treating A and B as constants 
and obtain the critical value of A. (The accuracy is extremely good because, in contrast 
to the determination of the UFB bounds, A and B are running very slowly at r c ). The 
value of A is dependent on m , but if we want to add enough .R-parity violation to remove 
the UFB problem entirely, we should do it at itlq = where the minimum is deepest. We 
therefore choose mo = in our evaluation of r c and also (for simplicity) Aq = 0. We find 
critical values 

0.76 > \ crit > 0.50 (78) 

with the upper value corresponding to tan/3 = 1.8 and the lower value to tan/5 = 15. (It 
drops to a plateau in tan/3 in much the same way as figure (4).) The revised estimate for 
the .R-parity violation required to remove all UFB bounds is therefore 

, 0.76h T M 1/2 Mi/a 
A > « 0.017^^ 

. 0.26h b M V2 Mi/2 , s 

A'> ^-^^ 0.009^^. 79 

\l (X 

Given that it is the a 6 terms which lift the flat directions, what set of -R-parity violating 
couplings do we need? For the two C directions we can get away with having only one large 
coupling, A123. This coupling can satisfy Eq. (|79"D without violating current experimental 
bounds assuming that fx is not very much smaller than M1/2] charged current universality 
requires that 

^123 < 0.05 f ) • (80) 



,100 GeVJ 

For the K directions, we need at least three relatively large A' couplings. In addition if we 
were, for example, to choose them to be A' 113 , A 213 and A 313 , then it is clear that we could 
just rotate the leptons into a basis where this is equivalent to just one coupling and we 
would still be left with two F-flat directions. We need to choose either A' la3 , \' 2b3 , A' 3c3 7^ 
or A' 13a , A' 23b , A' 33c 7^ 0, where a, b, c are all different (or any choice of three of these where the 
last two indices are all different). As an example, we choose A' 113 , A 223 , A' 333 7^ 0. The last 
two of these couplings are hardly constrained at all, and can easily satisfy experimental 
bounds and Eq.([73|) simultaneously. The first one however appears in quite restrictive 



bounds involving products of itself with one other A' coupling. It is also constrained by 
charged current universality to be 

x ^ <om (wtv)- (81 > 

Applying Eq . (|79|) to the product bounds in Ref. |15| results in new bounds for some of the 
other .R-parity couplings. Assuming that the squark masses are 100 GeV, for consistency 
with Ref. [O , they are 
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coupling <, 


Urrs limit y^\[i/ivii 2) 


• 

experimental limit 




U.UUUo 


U.UUUoo 


A 112 


D (14 




A 121 


0.03 


0.035 


A 122 


0.04 


0.02 


A 131 


9.0 x 10" 6 


0.035 


\/ 

A 132 


0.17 


0.33 


\/ 

A 213 


6.0 x 10" 6 


0.09 



When /i « Mx/2 we see that most of the new bounds are comparable to the current 
experimental ones. The bounds on A' 132 and A' 121 are actually more restrictive and those 
on A' 131 and A' 213 much more restrictive. 

To summarize, in minimal supersymmetry with constrained (degenerate) supersym- 
metry breaking, the addition of i?-parity violation just below the current experimental 
bounds (c.f. Eq.(|T9"D) is enough to lift dangerous flat directions and destroy local minima. 
For values of tan f3 which are away well from the quasi-fixed point the magnitude of the 
imparity violating couplings required to do this drops rapidly by a factor of ~ 2. 

5 Conclusions 

Using approximate analytic solutions to the RGEs of the MSSM, we have examined the 
formation of minima in F and D flat directions. We have given arguments which suggest 
that the traditional UFB bounds are not meaningful since, even if a radiatively induced 
minimum is not global, the decay rate back to the physical vacuum is very small. We have 
rederived the traditional UFB bounds analytically, and shown that they are numerically 
very close to satisfying a more meaningful sufficient condition - that there be no radiatively 
induced minima along the flat direction at all. The sufficient condition was derived for 
all tan/? < 15. 

In the literature the question of dangerous flat direction is usually tackled by either 
making an implicit or at least not very detailed appeal to cosmology, or by completely 
excluding regions of parameter space where the physical vacuum is metastable - a con- 
dition which, as we have seen, is probably not even relevant. The unbiased approach we 
have been advocating is to locate in the parameter space regions where some (unknown) 
cosmological input is needed to save the model. The 'traditional' UFB bound does not do 
this since it is still possible to get caught in a charge and colour breaking minimum even 
if it not the global one (and hence satisfies the 'traditional' UFB bound). The sufficient 
condition we gave is only slightly more restrictive than the 'traditional' bound but is in 
our opinion the only possible relevant condition and is very different in spirit. 

We also examined the effect of i?-parity violation, and found that a small amount of 
i?-parity violation (consistent with but just below current experimental limits) is able to 
lift the most dangerous flat directions and remove radiatively induced minima. In this 
way i?-parity removes the metastability problems which plague the MSSM, and rescues 
some models which could be worth saving. One example is the MSSM at the quasi-fixed 
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point which has much less dependence on initial parameters (as we have seen) and hence 
is more predictive |13 , |6|. This model is likely to be severely constrained by a combination 
of UFB bounds and higgs bounds coming from LEP2 [13[ |16| . A second example is the 



dilaton breaking scenario which has degenaracy in the supersymetry breaking terms (at 
tree level) and hence naturally suppressed flavour changing neutral currents. This model 
was excluded because of UFB bounds in Ref. IH7I . 



6 Acknowledgements 

We would like to thank Ben Allanach and Sacha Davidson for discussions and for a critical 
reading of the manuscript and also Toby Falk for additional comments. 



22 



References 

[1] J.-M. Frere, D.R.T. Jones and S. Raby, Nucl. Phys. B222 (1983) 11; M. Claudson, 
L. Hall and I. Hinchcliffe, Nucl. Phys. B228 (1983) 501; H.-R Nilles, M. Srednicki 
and D. Wyler, Phys. Lett. B120 (1983) 346; J-R Derendinger and C. A. Savoy, 
Nucl. Phys. B237 (1984) 307; R Langacker and N. Polonsky, Phys. Rev. D50 (1994) 
2199; A. Kusenko, R Langacker and G. Segre, Phys. Rev. D54 (1996) 5824; H. Baer, 
M. Bhrlik and D. Castano, Phys. Rev. D54 (1996) 6944 

[2] H. Komatsu, Phys. Lett. B215 (1988) 323. 

[3] J. A. Casas and S. Dimopoulos, Phys. Lett. B387 (1996) 107; J. A. Casas, 
hep-ph/9707475 ; J. A. Casas, A. Lleyda and C. Munoz, Nucl. Phys. B471 (1996) 
3. 

[4] J. A. Casas, A. Lleyda and C. Munoz, Phys. Lett. B389 (1996) 305 

[5] T. Falk, K. A. Olive, L. Roszkowski and M. Srednicki, Phys. Lett. B367 (1996) 183; 
A. Riotto and E. Roulet, Phys. Lett. B377 (1996) 60; A. Strumia, Nucl. Phys. B482 
(1996) 24; T. Falk, K. A. Olive, L. Roszkowski, A. Singh and M. Srednicki, 
Phys. Lett. B396 (1997) 50 

[6] S.A. Abel and B.C. Allanach, Phys. Lett. B415 (1997) 371. 

[7] F. Buccella, J-P. Derendinger, C. A. Savoy and S. Ferrara, Phys. Lett. B115 (1982) 
375; R. Gatto and G. Sartori, Phys. Lett. B157 (1985) 389; M. A. Luty and 
W. Taylor, Phys. Rev. D53 (1996) 3399; T. Ghergetta, C. Kolda, S. P. Martin, 
Nucl. Phys. B468 (1996) 37. 

[8] A. Linde, phys. Lett. B70 (1977) 206; 100B (1981) 37; A. Guth and E. Weinberg, 
Phys. Rev. D23 (1981) 876. 

[9] S. Coleman, Phys. Rev. D15 (1977) 2929; S. Coleman and C. Callan, Phys. Rev. 
D16 (1977) 1762. 

[10] M. Dine, W. Fischler and D. Nemeschansky, Phys. Lett. B136 (1984) 169; O. Berto- 
lami and G. G. Ross, Phys. Lett. B183 (1987) 163; M. Gaillard, H. Murayama 
and K. A. Olive, Phys. Lett. B355 (1995) 71; G. Dvali, Phys. Lett. B355 (1995) 78; 
M. Dine, L. Randall and S. Thomas, Phys. Rev. Lett75 (1995) 398; Nucl. Phys. B458 
(1996) 291. 

[11] J. E. Ellis et al, Nucl. Phys. B373 (1992) 118; S. Sarkar, Rept. Prog. Phys. 59 (1996) 
1493. 

[12] C. T. Hill, Phys. Rev. D24 (1981) 691; C. T. Hill, C. N. Leung and S. Rao, 
Nucl. Phys. B262 (1985) 517 

[13] V. Barger, M.S. Berger, P. Ohmann and R. J. N. Phillips, Phys. Lett. B314 (1993) 
351; V. Barger, M.S. Berger and P. Ohmann Phys. Rev. D49 (1994) 4908; M. Carena, 
M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. B419 (1994) 213; 



23 



426 (1994) 269; M. Carena and C. E. M. Wagner, Nucl. Phys. B452 (1995) 45; 
M. Lanzagorta and G.G. Ross, Phys. Lett. B349 (1995) 31; M. Lanzagorta and 
G.G. Ross, Phys. Lett. B364 (1995) 163. 

[14] K. Inoue, A. Kakuto, H. Komatsu and S. Kakeshita, Prog. Theor. Phys. 68 (1982) 
927; Errata ibid. 70 (1983) 330. 

[15] H. K. Dreiner, hep-ph/9707435 ; G. Bhattacharyya and A. Raychaudhuri, 
hep-ph/9712245 ; G. Bhattacharyya, hep-ph/9709395 

[16] J. A. Casas, J. R. Espinosa and H. E. Haber, hep-ph/9801365 ; S. A. Abel and 
B. C. Allanach, hep-ph/9803476 

[17] J. A. Casas, A. Lleyda and C. Munoz, Phys. Lett. B380 (1996) 59 



24 




25 




Figure 2: The value of m§ corresponding to an extremum appearing at that field value 
(parameterized by r). Each value of thq generates one maximum (to the left of the cross) 
and one minimum (to the right). This example is at the quasi-fixed point. 
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Figure 3: The height of the extrema vs their field value (parameterized by r). This 
example is at the quasi-fixed point and corresponds to substituting the wig values in 
figure 2 into the potential. 
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m /M 1/2 




Figure 4: The sufficient condition for tan/3 < 15 and for any value of Aq. A minimum is 
formed radiatively when m is below the line. 
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